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I. INTRODUCTION 

There is a revival of interest in the light hadron spec- 
troscopy. Several experiments are currently carrying out 
searches for new and exotic resonances and more facilities 
are planed for the near future. This is because the current 
knowledge of the spectrum is still very poor; in particular 
in the meson sector the Particle Data Group [T] lists ap- 
proximately a dozen or so of well established resonances 
and a plethora of states with poorly determined charac- 
teristics whose existence, in some cases, is even question- 
able. Establishing resonance parameters requires partial 
wave analysis and in particular knowledge of amplitudes 
in an unphysical region of the energy. This necessitated 
use of methods, such as dispersion relations, that ex- 
plore analytical properties of the physical amplitudes. 
Since the unitarity condition relates the measured cross 
sections to the discontinuities of the amplitude, disper- 
sion relations can be solved to obtain amplitudes for any 
complex (unphysical) value of the energy. Unfortunately, 
even in the simplest cases, no complete information about 
discontinuities of the amplitude is available and approxi- 
mations must be devised to solve the dispersion relations. 
Nevertheless, application of dispersion relations can lead 
to precision studies of resonance parameters, as was re- 
cently shown to be the case of the a resonance in 5- wave 
TTTT scattering [2H7]. Not only that but, exploring the 
relation between chiral effective field theory and QCD 
parameters, it was possible to shed new light onto the res- 
onance structure [SHU]- Further insight into the micro- 
scopic nature of a resonance can be obtained by studding 
its shape in different production process [T3HT5] . For ex- 
ample, a spatially extended source that produces a pair of 
mesons is not expected to couple strongly to a resonance 
that is dominated by valence quarks, i.e. whose wave 
function is spatially compact. Thus a compact source 
such as charmonium, when decaying into light mesons, is 
expected to couple to short-range components of inter- 
mediate states and resonances. 

It thus follows that in partial wave analysis it is im- 



portant to use amplitudes with proper analytical behav- 
ior across both the physical (right) and unphysical (left) 
cut. The discontinuity of the amplitude across the right 
hand cut (r/ic) is constrained by unitarity and relates 
partial wave amplitudes in production and scattering re- 
actions. The latter "contains" resonances and these will 
also appear in the production process albeit with mod- 
ified characteristics, that depend on the nature of pro- 
duction dynamics. The left hand cut (Ihc) determines 
the production dynamics and may enhance or suppress 
a particular resonance. For example, the QCD nature of 
the fli as a member of the quark multiplet was established 
by application of analytical methods to the Stt produc- 
tion in np collisions |16| . More recently, an analysis of 
the TTTT production in 77 fusion based on dispersion rela- 
tions showed that the a meson is expected to contain a 
substantial short range component, besides the dominant 
TTTT component 17J. The long range component originates 
from long range interactions between pions and, only in 
the limit of asymptotically large number of colors, it is 
excepted to be suppressed compared to the quark com- 
ponent. 

Unfortunately it is a common practice to ignore the 
above mentioned intricacies of amplitude analysis and 
rely on simple parameterizations i.e. in terms of a su- 
perposition of Breit-Wigner (BW) resonances. More so- 
phisticated analyses would use pole parameterizations for 
the iiT-matrix and implement unitarity restrictions across 
the right hand cut, but the role of the left hand cut is 
often overlooked. The main motivation of this paper is 
illustrate its role using as an example the tttt P-wave 
amplitude which is dominated by the well-known p(770) 
resonance. We will show how the general analytical prop- 
erties of amplitudes discussed above emerge in this par- 
ticular case and how the behavior of the amplitude on 
the Ihc can modify the resonance production. 

The paper is organized as follows. In the next section 
we summarize the analytical properties of partial wave 
amplitudes, both in the case of scattering and produc- 
tion. In Section |III| we solve the dispersion relation for 



the P-wave amplitude and discuss a simple parameteri- 
zations showing how the shape of the p-meson changes 
depending on the characteristics of the production pro- 
cess. Summary and outlook are given in Section [TV] 



II. ANALYTICAL PROPERTIES OF 

SCATTERING AND PRODUCTION PARTIAL 

WAVE AMPLITUDES 

In this work we focus on partial wave amplitudes which 
are a function of a single energy variable. This applies, 
for example to photoproduction of pion pairs on the nu- 
cleon, 7p -> vr+Tr^p. At fixed photon energy Ey, momen- 
tum transfer, t between the target and recoil nucleon and 
photon and nucleon helicities A^, the di-pion production 
amplitude, ^ is a function of the di-pion invariant mass 
squared s, and the spherical angle il = {6, (p) which de- 
scribes the direction of motion of the tt"*" in the di-pion 
rest frame. 



A = A{s,n,E-^,t,X,). 



(1) 



The angular dependence can be expanded in a series of 
partial waves labeled by the spin of the di-pion pair, I and 
its projection on, e.g. the photon direction ( t-channel 
helicity frame) 



A = ^Ai,„(s;S^,i)i1™(f^). 



(2) 



Finally the charged di-pion state can be decomposed into 
states of total isospin, 7 = 0,1,2 
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(3) 



The partial wave analysis can now be performed in bins 
oi Ej and t for a set of partial wave amplitudes Aimjia of 
a singe variable s - the di-pion mass squared. The singu- 
larities of these amplitudes in the complex-s plane will in 
general depend on t. This parametrization was recently 
used in the analyses of the two pion photoproduction 
data from CLAS at JLab [TFl [TO] and is currently un- 
derway for the hadronic decays of light charmed hadrons 
using CLEO and BES data gUj- 

In the following we will concentrate on production of 
pion pairs. The case of heavier mesons and/or baryons 
can be formulated analogously, albeit with complications 
arising from presence of sub-threshold cuts and spin. As 
a function of the di-pion invariant mass squared, s ~ 
4((7^ -|- m?), with q being the relative momentum, and 
at fixed values of other kinematical variables, the l-th 
partial wave production amplitude is factorized into the 
angular momentum barrier factor, q^ and the amplitude 
F{s) that is free from kinematical singularities 



A{s) = q^Fis). 



(4) 



The amplitude F{s) is a real analytical function with a 
right hand cut starting at threshold, Sth = 4m^ and with 



each open inelastic channel, Si > Sth contributing to the 
discontinuity across the cut. The value of this disconti- 
nuity is constrained by unitarity. For s on the positive 
real axis, defining F^{s) = F{s ± ie) as the boundary 
value of the function F{s) on the upper (lower) lip of the 
cut, unitarity implies 

ImF+is) = t-is)p{s)F+{s)e{s-Sth)+CTis)eis~s,). (5) 

Here t~ (s) = i(s — ie) is the boundary value on the lower 
lip on the right hand cut of the tttt 1-wave scattering 
amplitude t{s) considered as an analytical function of s. 
The two body phase space is given hy p = (1 — Sth/s)^ 
and cr{s) represents the contribution from production of 
inelastic channels above Si. We are using the following 
normalization of the tttt amplitude. 



r){s)e^'^'^'^ = l + 2ip(s)t+(s), 



(6) 



where 5, 77 is the phase shift and inelasticity, respectively. 
In addition to the right hand unitarity cut, the produc- 
tion amplitude is discontinuous for negative s. The lo- 
cation of the left hand cut s < sq < depends on the 
underlying production dynamics. It is related to thresh- 
olds for particle production in the crossed channels and 
plays the role of the driving term in the integral equa- 
tion for the amplitude that follows from a dispersion re- 
lation. It has a similar role to that of the potential in 
the non-relativistic Lippmann-Schwinger equation of the 
Schrodinger theory. Thus for s real and s < sq. 



ImF+{s) = ImFL{s + ie)6'(so - s), 



(7) 



with the "potential" Fl{s) defined as a real analytic func- 
tion in the complex s-plane with the Ihc discontinuity 
given by Eq.Q, 



1 r° ImF+{sL) 

Fl(s) = - / dsL 

TT J-00 sl - s 



(8) 



From Eqs.([5]),([7| it follows that F{s) satisfies an integral 
equation. 



F(s) = FLis) + - I dsB.- 



■n 



Sth 



sr- s 



T Js, Sr- S 



(9) 



where we assumed that the integrands vanish as \s\ — > cx3. 
If F{s) is bound by a polynomial in s, the dispersive in- 
tegral can be made convergent by subtractions. These 
introduce additional parameters, that are related to the 
asymptotic behavior of the amplitude. In general little is 
known about the left hand cut discontinuity, i.e Fl{s). 
For example in bootstrap calculations it is approximated 
by particle, or more generally Regge exchanges with pa- 
rameters adjusted so that the solution of the dispersion 
relation reproduces the known resonances. More recently 



chiral effective field theory has been used to construct ap- 
proximation to Fl at low energies by expanding it in pow- 
ers of s/A^, where A^ — Anfj^ is the chiral scale [31I51[^H- 
[25] . The dispersion relation in Eq. ^ is an integral equa- 
tion for the production amplitude F(s) which as input 
takes i) the scattering amplitude t, ii) the amplitude 
representing contributions from production of inelastic 
channels a and Hi) the "potential", F^. The analytical 
solution of Eq.Q is known, and involves the function 
D{s) from the N/D decomposition of the scattering am- 
plitude, t{s) piH25] . Even though the N/D method for 
solving partial wave dispersion relations has been exten- 
sively studied in the past, due to its important role in 
partial wave analysis we summarize its main features in 
the following paragraphs. 



A. N/D representation of the scattering amplitude 

The dispersion relation for the partial wave scattering 
amplitude t(s) is given by 



SL 



SR 



m 

It becomes an integral equation for i(s), once unitarity is 
implemented. It relates the rhc discontinuity in t to the 
amplitude itself via. 



Imt+{sR) = p{sR)\t{sR)\'^0{sR- sth) 
, l-V^{sR.) ^, . 

+ -. &{SR-S,). 



(11) 



The solution of Eqs.( 10 1,(11 1 is expressed in the form. 



t{s) = (S- Sth) 



D{s)-^ 



(12) 



where the angular momentum barrier factor has been ex- 
plicitly factored out and N{s), and D{s) are defined so 
to have only left and right hand cut, respectively. The 
definition of N and D is unique up to an overall constant, 
which we fix by normalizing D at the elastic threshold, 
D{sth) = 1. For / > 1 the uniqueness follows from the 
threshold behavior, t{sth) = and the asymptotic behav- 
ior at large s, \t{oo + ie)\ < 0(1). For S-waves, one sub- 



traction may be needed to make the integrals in Eq.(lO) 



convergent and consequently the dispersion relation for 
t (or N or D) contains one undetermined constant {i.e. 
scattering length). 

From now on we focus on the P-wave, / = 1 amplitude. 
From the dispersion relation for t it follows that. 



7V(s) = 



1 



'" , Imt(sL)D(sL) 

asL-, r-, r 

-oo {sl~ Sth){SL- s) 



D{s) = 1 



Sth 



ds, 



p{sr)N{sr)R[sr) 



SR 



(13) 



where R = /mt"'"/p|ip, differs from unity for sr > Si 
because of inelastic channels opening. Substitution of 
the equation for A'' into the equation for D leads to an 
integral equation. 



Dis) 



where the kernel 



dsLK{s,SL)DisL) (14) 



K{s,si 



Sth 



dsR 



p{sR)R{sR)Imt{sL) 

{sL - Sth)iSL - Sr){s]i - s) ' 

(15) 
is given in terms of the "potential", i.e discontinuity of 
t across the left hand cut, and the inelastic contribution 
proportional to i?— 1. Since |i(cx)-|-ie)| < 0(1), it follows 



from Eq.(13) that asymptotically |A^(s)| < 0{l/s) and 
|I?(s)| < 0(1) and no subtractions are needed. 

The dispersion relation for D in Eq.([l3| assumes i 7^ 0, 
which is always the case in the inelastic region when 
?7 7^ 1. Otherwise zeros of t correspond to poles of D, 
the so called CDD poles that are not accounted for by 
Eq.(13) [27]. Therefore in presence of zeros of the scat- 
tering amplitude, the CDD poles have to be added "by 



hand" and result in, 



Dis) = l 



S- Sth 



N„ 



-(S - Sth)IlpZl 



ds, 



7p 



P{sr)N{.Sr)R{.Sr) 



SR- S 



^ [Sp - s){Sp ~ SthY 



(16) 



where 7^ and Sp are the (real) residue and position of 
the p-th CDD pole, respectively. At every CDD pole the 
phase of the amplitude passes through 180° and Levinson 
theorem relates the number of CDD poles to the phase 
shift at infinity, ^(oo)/7r = Np [55] . If the residue of a 
CDD pole is small, then D{s) will develop a zero on the 
second sheet near the position of the pole, i.e. will pro- 
duce a resonance. Thus, in the past it has been proposed 
to identify CDD poles with the elementary quark bound 
states that turn into physical resonances when coupled to 
the continuum channels. Indeed it has been shown that 
in potential models describing, for example scattering off 
a static source with internal structure, CDD poles in the 
dispersion relation for a scattering amplitude correspond 
to excitations of the target. There is, however no proof of 
such correspondence in QCD, and as we will show in the 
following section, a correspondence between CDD poles 
and quark model states may be more complicated than 
in potential theory. 

It follows that, even though unitarity and the left hand 
cut discontinuity do not yield a unique solution of the dis- 
persion relation for the scattering amplitude i, the CDD 
ambiguity disappears if the phase shift and inelasticity 
are known. In this case, the solution of Eq.(16) is given 
by, 



D{s) = n 



p=l 



Sth 



^<t>{s), 



(17) 



where the first factor comes from the ODD poles and the 
Omnes-Muskhehshvih function il 1291 130] 



O^ = exp - 



S - Sth 



(sr) 



dsR-. — ,, 
[sr-s){sr 



Sth), 

.(18) 

is given in terms of the phase of t, i = |t| exp(i0), which is 
identical to the phase shift 5 in the elastic region. Outside 
the elastic region the phase (j) is determined from phase 
shift and inelasticity and is bound between and tt. 



B. N/D representation of the production 
amplitude 

The solution of the dispersion relation, Eq.Q for the 
production amplitude can be represented in a number 
of equivalent ways f31| . From the point of view of data 
parametrization, a particularly useful representation is 
given by 

F{s) = ^ = F{stH) + GL{s)-G.{s) ^ ^^g^ 



D{s) 



D{s) 



where D{s) is given by Eq.(17) (or (16l), the numerator 
function Gl contains only the left hand cut, 



Gl{s) = 



Sth 



ImFL{sL)D{sL) 



dsL-, \.^ ^ ', (20) 

■ oo {SL~ S)[SL~ Sth) 



while the function Gi{s) describes inelastic contribution 
to production, 



where s — > zi^(^r){s) represents a conformal mapping of a 
cut plane on a unit circle |32j . 



zl{s) 
zr{s) 



1 - y/l - s/Sl 
1 + v/1 - S/SL 
1 - y/l - S/S, 

1 + v/1 - s/s, ■ 



(23) 



One would then determine the coefficients Ci^L,Ci,R from 
fitting F{s) to the data. 

Another way of representing the solution of Eq.(|9|, 
which in particular amplifies the spacial characteristics 
of the production process, is obtained by represent- 
ing D{sl) in Gl{s) through its dispersion relation (c./. 
Eq.(16)), resulting in 



F{s) = Fl{s) + 



F{sth) - FListh) 



Sth 



ttD{s) 

S - Sth 

TrD{s) 



dsR 

[SR 

) 

dsj 



D{s) 
ImD{sR)FL{sR) 



s){sr- Sth) 

Ima{sR)D{sR) 



t {sr)p{sr){sr- S){SR- Sth 



(24) 



Both representation, Eq.(19l and Eq.(24), below the 



inelastic threshold s < si, satisfy the final state inter- 
action theorem, arg F{s) — — arg D{s) = arg t{s). In 
Eq. ( 24 ) , the production amplitude is expressed in terms 



of the production amplitude Fl evaluated in the physical 
{s > Sth) kinematics. Parametrizing F^ via the confor- 
mal map {c.f. Eq.(23)), 



FUs)=J2c,zUs), 



(25) 



G,(s) 



Sth 



dsj 



Ima{sR)D{sR) 



"t {sr)p{sr){sr- S){SR- Sth)' 

(21) 
We have assumed that at most one subtraction is needed, 
which is the case if D has a CDD pole at infinity that 
as we will discuss later is the likely case for the P-wave. 
Expression for the production amplitude in Eq.(19l can 
easily be adopted to fits of experimental data. Assum- 
ing that D{s) is a known function (we will discuss spe- 
cial cases in the next Section), the numerator in Eq.(19) 
depends on the inelastic contribution (cr) and the pro- 
duction dynamics {Fl). These are the quantities one 
should extract from the data. The form in Eq.(19) sug- 



gests that, instead of parametrizing Fj^ and a, it would 
be more efficient to parametrize the entire expression ap- 
pearing in the numerator. It represents a sum of two 
analytical functions one with the left an the other with a 
right hand cut that begins at the first inelastic threshold. 
Thus one can conveniently use a series expansion. 



Gl{s) 



i=0 



Ct.LZLis), Gi{s) 



i=0 



c,^Rz}^is), (22) 



and replacing the inelastic contribution with the fit func- 



tion Gi from Eq.(23) leads to 



F{s) = 



F{sth) 



4(s) 



D{s) '^ 
with the functions Ii{s) given by 



z'jsth) 
D{s) 



his) 
D{s) 



G^js) 

' Dis) 

(26) 



his) = - 



S - Sth 



ImD{sR)zl{sR) 
dSR- r. (27 j 

[SR- S)[Sr- Sth) 



The asymptotic behavior is encoded in the rate of conver- 
gence of ^ , CjZf at the circle \zi^\ = 1. The applicability 
of Eq. ( |26| , in experimental data fit is however based on 
the assumption that only a few orders in the expansion 
in powers of zj^ are needed to describe a particular data 
set, i.e one wants the fits to be insensitive to the poorly 
known asymptotic behavior. 



In the elastic region sth < s < Si from Eq.(24) one 
finds 



F{s) 



^ giH'^) 



Fl(s)cos(5(s)-S(s) 



sin^(s) 
^(^ 



(28) 



where S is the elastic phase shift and i? is a real function 
given by 



B{s) 



FL{sth) ~ F{sth) 



t{s)D{s) 
s - sth p y [°° , ImD{sB)FL{sR) 

■Kt{s)D{s) ■ ' Jg^^ ^{SR- S){SR- Sth) 



S - Sth 



dsR- 



Ima{sR)D{sR) 



Trt{s)D{s) J^. t {sr)p{sr){sr - s){sr - sth) ' 

(29) 



with P.V. standing for the principal value. In particular, 
if at some s — Sr the phase shift passes through a res- 
onance, i.e QO'^, the first term in Eq.(28l will produce a 



zero in the amplitude, while the second term will produce 
a peak there. The zero is due to a destructive interfer- 
ence between the direct production of two-particles and 
their final state interaction. In the case were production 
is short-ranged, as discussed earlier, the start of the left 
hand cut in Fl{s), sq is far away from the physical re- 
gion, S{) « Sth, the principal value integral extends over 
a large interval in sr, sr < |so| and the contribution from 
the B term in Eq.(28) is significant. On the other hand. 



if the production source is diffuse and inelasticity is small 
the resonance peak may be significantly distorted by the 
zero from the cos S term. Thus analyzing the shape of 
a resonance can shed light on its production character- 
istics and thus its nature. We also note that the first 



term in the rhs. of Eq.(29) contributes if the dispersion 
relation for G{s) requires subtraction. It is needed if 
D{s — > oo) ^ 0{s) i.e D{s) has a CDD pole at infinity. 



III. P-WAVE TTTv SCATTERING AND 
PRODUCTION 



The solution of the dispersion relation for the scat- 
tering amplitude requires knowledge of the inelasticity 
{R — 1) and the left hand discontinuity. In the case of 
the P wave the former is know up to y/s = 1.9 GeV and 
the latter has been evaluated in Ref. [33]. In the case 
of TTTT scattering, crossing symmetry gives an additional 
constraint between the "potential" and the right hand 
discontinuity. In Ref. |33) a particular parametrization, 
for the left hand discontinuity 



Imt{s) 



bs 



+ c[l — cos(27ra;)] 



(30) 



with a = 0.48 GeV*, b = 1.21 GeV^ 



c = 0.601 and 
2^^o.36 f^^ , / _qo,^2 



X = (-0.320/(0.283 -I- s/ GeV^))'J■^^ for s < -32m^ has 
been shown to faithfully represent this constraint. For 



-32?Ti^ < s < So = 0, crossing symmetry leads to 



27772 rsth-s / 2sd 

S - Sth Jstj, V S- Sth 



x^Ci7^(2; + l)/TOiu(sfl)P/ (1 



/=o 1=0 



2s 



Sr - Sth 



(31) 



Since the integration range is limited between threshold 
and 36to^ ~ 0.675 GeV^, the sum can be truncated to in- 
clude only a few low partial waves and the corresponding 
amplitudes Imtij (t = iii) expressed in terms of known 
phase shifts and inelasticities. For R in Eq.(13l we use 
the data from |3il - l3S] and assume 7/ = 1 for ^/s> 2 GeV. 
The kernel in Eq.(15) is then inverted numerically and 



the D function is obtained from computing 



D{s) 



dsL[l + K]-\s,SL)P{sL) (32) 



where P{sl) ~ 1 if no CDD poles are assumed, or 
P{sl) = l-l-7p(sL-Sfh)/(sp-St,i)(s-sp) forthecaseof 
one CDD pole. Since the P-wave phase shift stays below 
180° all the way up to the first relevant inelastic thresh- 
old associated with KK production, no CDD poles are 
expected. The exception is a pole at infinity. In our our 
case, since we have assumed that inelasticity is negligible 
at large s, the CDD pole at infinity could in fact appear 
at a finite sp > (2 GeV)^. Once D{s) for s < is ob- 



tained from Eq.(32), both N{s) and D{s) in the physical 



region can be computed from Eq.( 13 ) and compared with 



the measured phase shift and inelasticity. The results are 
shown in Fig. [T] respectively with the dashed-dotted line 
showing the result computed without the CDD pole and 
the solid with one CDD pole, whose residue and pole po- 
sition have been fitted to the scattering data [34ll36] . It 
is clear that the resonance behavior corresponding to the 
p(770) is due to the presence of a CDD pole. The exact 
location of the pole cannot be established due to the un- 
known behavior of the phase shift at high energies. In the 
current fit we find the CDD pole with 7p = 50.4 GeV^ 
and Sp = 6.9 GeV^ which for all practical purposes as far 
as the p meson is concerned might be taken to be at in- 
finity. We also note that in |3] it was erroneously stated 
that the P-wave has a CDD pole at threshold. Vanishing 
of the amplitude at threshold is due to the angular mo- 
mentum barrier and not to the presence of a CDD pole, 
the latter being related to the dynamics of the scattering 
process. 

Next explore the connection between the CDD pole de- 
scription of the p resonance with that of the quark model. 
From the quark model perspective, the p meson is consid- 
ered as a quark-antiquark bound state which becomes a 
resonance when coupling to the continuum tttt channel is 
allowed. Generically, with such dynamical assumptions, 
the P-wave tttt scattering amplitude can be expressed as 
a solution of a separable integral equation shown in Fig. [2] 



a integral given by 




%(«) = - / dsR ^ ^ . (35) 

TT Js,^ SR- S 



In therms of the N/D representation (c./Eq.(12|) the 
ampHtude of Eq.(33l is becomes, 



Dis) = 1 



S - Sth 



' p{sr)N{sr) 

dsR (s - SthjC, 

sr- s 



N{s) = Cfis) (36) 

where the (positive) constant C is defined by 

C^^=m%-sth / dsRp{sR)f{sR) (37) 



Comparing with Eq.(16) shows that the quark model 




IGeVJ 



FIG. 1: P-wave phase shift (upper panel) and inelasticity 
(lower panel). Data from |34H36| . dahsed-dotted (solid) line 
solution of dispersion relation without (with) a CDD pole. 
Dashed line is the fit of the quark model from Eq.(33l. 



which leads to 



t{s) 






(33) 



Here f{s) is a vertex function which represents the cou- 
pling of the bare, quark model p state to the tttt channel, 
i.e. is given in terms of the overlap of the quark model 
p meson and the pion wave functions. In general /(s) is 
analytical in the s-plane, except for "potential" -like sin- 
gularities for negative s. In the pole approximation it is 
then given by (s/ > 0), 



/(^) 



XfS 






(34) 



The contribution from the open tttt channel results in the 
modification of the bare p propagator {m^ — s)~^ by loop 







FIG. 2: Quark model representation of the p resonance. The 
bare state with unrenormalised mass tub couples to a single, 
two-particle, open channel. 

representation of the p meson, in the absence of higher 
excited (p') states corresponds to a single CDD pole at 
infinity i.e. with Sp,7p — ?► oo, Jp/s^ = fixed. We fit 
the three constants Sf , C, Xf to the P-wave phase shift 
below -^s = 1.2 GeV, i.e. below the energy at which in- 
elasticity becomes sizable (since t in Eq. ( 33 ) describes a 
single channel). The result of the fit is shown in Fig. [I] 
by the dashed line. Below ^/s = 1.5 GeV the solution 
of the dispersion relation and the simple quark model 
parametrization are essentially indistinguishable. The 
comparison between N, D and t from dispersion relation 
calculation and the quark model is shown in Fig. [31 |4j 
and [5J respectively. For s > and in the p region, t 
computed from dispersion relations and the model should 
agree since both reproduce the phase shift and inelastic- 
ity. The D function is sensitive to the phase of the ampli- 
tude over the whole energy range and since the phases of t 
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FIG. 3: Imaginary (upper plot) and real (lower plot) of 
A^. Solid line is the solution of dispersion relation, Eq. (131 
, dashed line of the quark model, Eq.(33l. 



FIG. 4: Imaginary (upper plot) and real (lower plot ) of 
D"^ . Solid line is the solution of dispersion relation, Eq.( 13 1 



, dashed line of the quark model Eq.(33 1. 



A. p meson production 



in both cases approach 180" there is also little difference 
below the inelastic region. The only noticeable differ- 
ence is in the "potential", i.e. the numerator function 
N . In the case of the model in Eq.(33) there is a double 
pole at s = —Sf — —3.1 GeV" imposed by the simple 
parametrization of Eq.(|34|, while the actual solution of 



the dispersion relation shows a structure at s ~ — 1 GeV 
and s ^ 2.5 GeV^. It is clear that knowledge of the am- 
plitude in the physical region alone does not allow one 
to reconstruct the potential. This is because the p is 
dominated by the CDD pole and the "potential" inter- 
action, i.e. rescattering of the pions, has little effect on 
the p. This can also be seen in Fig. [l] where the compar- 
ison between the full solution and the no-CDD solution 
is shown. 



We finally illustrate the role of the left hand cut 
on the p shape in a specific production process. In 
a quark model description the production amplitude 

A{s) = qF{s) = (s - SthY^^Fis), {c.J. Eq.iH) becomes 



Fis)^9{s) + ^ 



fis)Igf{s) 



Iffi^y 



where 



hfi^) = - 



ds, 



{sR - sth)pisR)g{sR)fisR) 



sr- s 



(38) 



(39) 



and is depicted in Fig. [6] The first term in Eq. ( 38 ) cor- 
responds to direct production of the two-pions, and is 
given by a "potential" term, g{s), which we parametrize 







FIG. 6: Quark model representation of the p production 
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FIG. 5: Imaginary (upper plot) and real (lower plot) of 
t. Solid line is the solution of dispersion relation, Eq. (131 



, dashed line of the quark model, Eq.(33l 



5(5) 



^3^3 



(40) 



with Sg > characterizing the range of interaction or the 
size (^ l/\/Sg) °f ^^"^ source. A compact, or point-like 
source corresponds to g(s) — >■ const., i.e. Sg — >■ oo while 
a spatially extended source has s^ — >■ (with XgSg fixed). 
The production amplitude in Eq.(|38|) can be cast in the 



form of the general representation of Eq. ( 19 1 . In this case 
Gi{s) = since the model contains no inelastic channels, 
and 

G{s) =. C [9{s){ml -s)- {9{s)Iff{s) - I{s)Igf{s))\ . 

(41) 
Indeed, the right hand cut singularities of /// and Igf 
cancel out and G{s) has only the left hand cut singular- 
ities of g{s) and /(s). Similarly Eq.(38) can be written 



in the form of Eq. ( 28 1 with 



/toFl(s) = ~T:5{s + Sg)\gSgFg-T:6{s+Sf)\fS}Ff, (42) 
where 

,m\ + Sg-Iff{-Sg) 



Fg ^ C 



Ff ^ C 



-•f 



Di^Sg) 



and so 



FLi.s)^Fgg{.s) + Fffis). 



(43) 



(44) 



In the case of a diffuse source the interaction range is 
small, and the principal value integral in Eq. ( 29 1 is sup- 



pressed. In this case the production amplitude is sup- 
pressed at the resonance mass, since the other term is 
proportional to cos S which vanishes. In the quark model 
model case, however, we see that regardless of the range 
of TTTT production, determined by g{s), the production 
potential also contains a term responsible for direct pro- 
duction of the quark state, i.e. the term proportional 
to /(s) in Eq.(44|. Thus the resonance associated with 



a quark state or CDD pole is expected to show up as a 
bump in a production process. As expected, the g term 
in Eq.(44) becomes more relevant as for diffuse sources, 
i.e. when s„ << st and P.V.I„f < P.V.Itf. This effect 



'•7 



is illustrated in Fig. [7j 



hf 



^ff- 



IV. SUMMARY AND OUTLOOK 

Development of analytical parameterizations of hadron 
production amplitudes is necessary for a successful par- 
tial wave analysis. Here we presented the analysis of 
the TTTT P-wave amplitude. The dominant feature in the 
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FIG. 7: Modulus squared of the production amplitude j-F(s)| 
from Eq.(38l computed for different values of the direct, tttt 
production range Sg. For large value of Sg, i.e Sg > 1 GeV^, 
the interaction range is small and the shape of a quark model 
(compact) resonance is unaffected. As the interaction range 
increases the rate of low mass pion production increases dis- 
torting resonance shape at lower energies. 



energy range, however for practical applications, it may 
be assumed that the CDD pole is at infinite energy, i.e. 
77(00) — > 1, (5(oo) — > TT. While this is consistent with the 
expectation for the asymptotic behavior of the pion form 
factor [33 [3S] it contradicts the expectation that inelas- 
tic open channels dominate at high energies (77(00) — )• 0) 
and requires further studies. The CDD pole in a par- 
tial wave amplitude is equivalent to presence of a " bare" 
state e.g. quark bound state, which couples to the tttt 
continuum. It is thus not surprising that a generic quark 
model description of P-wave tttt scattering with a single 
bound state faithfully describes the P-wave amplitude 
up to -^/s ^ 1.6 GeV where inelasticity becomes signifi- 
cant indicating a large KK component of the /9'(1600). 
While the p resonance shape can be distorted if the pro- 
duction amplitude varies strongly near threshold, to our 
best knowledge such features has not been observed in 
the experimental data. Even though the left hand cut 
behavior of the P-wave amplitude in scattering and pro- 
duction does not seem to drive the p resonance, this does 
not need to be the case for other resonances {e.g. the 
sigma meson) leading to the conclusion that the left hand 
cut properties of amplitudes require attention in partial 
wave analyses. 
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